#set pond parameters

	#must declare dimensions!
	pond=array(dim=c(2,5))	

	#column name #column number #description
	p.area=1 #area (ft^2)
	p.bottom=2 #bottom elev above MSL (ft)
	p.surf=3 #initial head (ft)
	p.min.surf=4 #min head (ft)
	p.max.surf=5 #max head (ft)
	
	pond[1, ]=c(area=1000,  bottom=1.97, surf=0,min.surf=0,max.surf=10) #Ox1
	pond[2, ]=c(area=1000,   bottom=1.31, surf=10,min.surf=0,max.surf=10) #Ox2

#set box parameters
	
	#must declare dimensions!
	box=array(dim=c(1,16))
	
	#column name #column number #description
	box.type=      	1 	#1=weir/sluice box #2=pump
	box.upstream=  	2 	#upstream pond
	box.downstream=	3 	#downstream pond
	box.crest=     	4 	#elevation of weir crest
	box.top=       	5 	#elevation of box top;
	box.width=     	6 	#width of weir crest
	box.roughness= 	7 	#empirical roughness factor (1=min, 0=max)
	box.pump=      	8 	#pumping rate
	box.l.up=      	9 	 #upstream pipe length
	box.d.up=      	10	 #downstream pipe length
	box.f.up=      	11	 #upstream pipe friction factor
	box.l.dn=      	12	 #upstream pipe diameter
	box.d.dn=      	13	 #downstream pipe diameter
	box.f.dn=      	14	 #downstream pipe friction factor
	box.gate.up=   	15	 #percent opening of gated pipe (upstream)
	box.gate.dn=   	16	 #percent opening of gated pipe (downstream)
	
	box[1,]=c(	#name='Ox1 to Ox2',
			type=1,
			upstream=1,
			downstream=2,
			crest=2,
			top=10,
			width=10,
			roughness=1,
			pump=0,
			l.up=10,
			d.up=2,
			f.up=0,
			l.dn=10,
			d.dn=2,
			f.dn=0,
			gate.up=0,
			gate.dn=1)
